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Entangled photons, generated by spontaneous parametric down-conversion 



^ . and image-processing applications. Since this source is an example of a three- 

. wave mixing process, there is more flexibility in the choices of illumination 



and detection wavelengths and in the placement of object(s) to be imaged. 
Moreover, this source is entangled, a fact that allows for imaging configura- 
tions and capabilities that cannot be achieved using classical sources of light. 



. In this paper we examine a number of imaging and image-processing configu- 

^ : 

^ . rations that can be realized using this source. The formalism that we utilize 



facilitates the determination of the dependence of imaging resolution on the 
physical parameters of the optical arrangement. 

I. INTRODUCTION 

The process of optical parametric three- wave mixing in a second-order nonlinear medium 
p|^IO| involves the coherent interaction between three optical fields with, generally, different 
wavelengths: pump, signal, and idler. Because of the phase-matching requirements [^|^,|TU[, 
the wave-vectors are related and the spatial distributions of the fields are therefore highly 
coupled. This process may therefore be utilized in distributed multi-wavelength imaging or 
image processing, where objects are placed in the path of one or more of these fields of 
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different wavelengths and the intensities or cross-correlations are measured. Two examples 
of phenomena based on three-wave mixing, optical parametric oscillation (OPO) and op- 
tical parametric amplification (OPA), have been studied extensively and many interesting 
phenomena of spatial correlation [|IT|, pattern formation |T2|,|13[, and reduced- noise image 
amplification [T^JT^ have been reported. 

A third example is the process of spontaneous parametric down-conversion (SPDC) [jl6|, 
a phenomenon that exhibits quantum entanglement fl^ . The signal and idler waves, created 
when the nonlinear medium is illuminated by an intense laser beam (pump), are produced 
in the form of photon pairs in an entangled quantum state (biphotons). Spatial and spectral 
entanglement are a consequence of the multiple possibilities for satisfying conservation of 
energy and momentum for each photon pair. Interest in the SPDC process has spurred 
many studies of its spatial and spatio-temporal photon correlation properties [[T8|-p^, and 
some imaging applications based on the measurement of photon-pair coincidence have been 
proposed [|29|,|g and tested |31|j3g. 

In this paper we develop a general Fourier-optics theory of image formation based on the 
SPDC process. In Section 2 we explore new configurations for multi-wavelength distributed 
imaging and image-processing applications. We follow an approach introduced in a previous 
paper in which we established a duality between partial quantum entanglement and 



classical partial coherence theory. We use the formalism developed in Ref. [ |3^] and apply it 
to the image-formation process in distributed multi-wavelength imaging configurations made 
possible by the nature of this SPDC source. In Section 3 we study the imaging resolution 
of these configurations and the effect on it of the various physical parameters of the system. 
In the Appendix we provide a brief review of classical imaging theory in the framework of 
the optical bilinear transformation. 
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II. CONFIGURATIONS FOR SPDC BIPHOTON IMAGING AND IMAGE 

PROCESSING 



The principal function of an optical imaging system is to transfer the spatial distribution 
of some physical property of an object (transmittance, reflectance, or absorbance), via an 
optical wave, to a remote location where it is measured with a photodetector An 
image-processing system transforms one image into another with enhanced features, or ob- 
tains a new image from more than one image, such as the correlation of two images P3|-P7|. 
We examine here various configurations for imaging and image processing based on bipho- 
ton beams generated by SPDC. As mentioned above, the existence of three optical beams 
(pump, signal, and idler) allows us to construct imaging configurations that are not achiev- 
able with other single-beam optical sources. We may place an object that is to be imaged in 
either of the three beams. We may alternatively place more than one object in these beams 
and obtain the correlation of their images. 

All these configurations include two features. The first is that they are examples of 



distributed imaging p9| , |30| , |38| . In analogy to distributed computing, where the computation 
resources are distributed over a network, distributed imaging allows us to reallocate the 
imaging components from the particular path connecting the source to the object to be 
imaged. The second is that they allow the possibility of multi-wavelength imaging: the 
object may be illuminated with one wavelength whereas detection takes place at another. 
The use of other nonclassical sources of light in imaging has been shown to lead to 



noise reduction ||1^JT5|J39|J40|] . In this paper we direct our attention to various imaging 
configurations. The quantum nature of the SPDC source offers the additional advantage of 
noise reduction, but this is immaterial to the task at hand. 

In the process of SPDC, an intense laser beam (pump) illuminates a nonlinear crystal 
(NLC) with quadratic nonlinear susceptibility ||7|JT6||. Some of the pump photons disintegrate 
into pairs of photons (known traditionally as signal and idler), which conserve the energy 
and momentum of the parent pump photon. Consider the situation depicted in Fig. 11]. The 
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pump beam illuminates the NLC and the signal and idler beams are measured by the single- 
photon detectors Di and D2, respectively. We assume throughout a planar source and a 
one-dimensional geometry in the transverse plane for the sake of simplicity but without loss 
of generality. Optical systems, containing objects to be imaged and any optical components, 
may be placed in any of the three available beams. 

The signal and idler photons can be emitted from the NLC in a variety of configurations. 
The two photons may be emitted in two different and distinct directions, in which case each 
photon will pass through a different (and possibly remote) optical system; this configuration 
is denoted non-coUinear. The two photons may be emitted in the same spatial wave packet, 
the collinear case, but have some distinguishing characteristic, such as frequency or polar- 
ization, whereupon they effectively pass through different optical systems if the components 
are dispersive or polarization dependent. The two photons are detected in the same output 
plane in this case. Finally, the two photons may be emitted in the same spatial wave packet 
and have no distinguishing characteristic, the collinear degenerate case, and thus they pass 
through the same optical system and are detected in the same output plane. 

The coincidence rate of photon pairs at the two detectors, Di and D2 located at positions 
Xi and X2, respectively, is proportional to the fourth-order correlation function of the fields, 
G^'^\xi,X2) p3i^ , the biphoton rate. The signal and idler beams traverse optical systems 
described by their impulse response functions hs{xi,x) and hi{x2,x), respectively. It has 
been shown that the biphoton rate is given by |Q 



G'(2)(xi,X2) = |V(Xi,X2)P, (1) 

where the biphoton amplitude ■ip{xi,X2) is 



here Ep{x) is the spatial distribution of the pump field at the entrance to the NLC. The 
result in Eq. (j^) was derived assuming a thin NLC and the presence of narrowband spectral 
filters in the optical system. These two assumptions simplify the analysis considerably 



without overshadowing the physics of the imaging processes discussed. They will be relaxed 
in Section 3. 

An interpretation of Eq. (0) that is useful in understanding the behavior of such a system 
was advocated by Klyshko [0,^,^,^ under the name " advanced wave interpretation" . In 
this picture, the biphoton amplitude in Eq. can be viewed as the impulse response 
function of an optical system represented by the cascade of three systems: propagation 
from Di at xi back through a system with impulse response function hl{x,xi) = hs{xi,x), 
modulation by the pump field Ep, and subsequent transmission through a system with 
impulse response function hi{x2,x). An intuitive advantage can be gained by unfolding the 
system in this way as will become clear shortly. 

Two special correlation functions deriving from G^'^\xi,X2) in Eq. (|l|) are of interest: 
the marginal coincidence rate J*^^)(x2), and the conditional coincidence rate Iq'\x2), defined 
by 

I^^\x2) = Jdx,G^^\xuX2), (3) 

II,'\X2)=G^'\0,X2). (4) 

The marginal coincidence rate I^'^\x2) is proportional to the probability of detecting a 
photon at X2 by D2 when detector Di detects a photon at any location —00 < a^i < 00. The 
conditional coincidence rate l'^\x2) is proportional to the probability of detecting a photon 
at X2 by D2 when Di detects a photon at Xi = 0. 

We proceed to examine the five configurations that are possible with this optical source 
and explore their imaging and image-processing potential. 

A. Object in the signal (or idler) beam 

The generalized biphoton optical system described by Eq. (Q) permits the object to be 
placed in either the signal or the idler beams such that its transmittance (or reflectance) 
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modifies either of tfie impulse response functions hs or hi. Without loss of generality, we 
assume that the object is placed in the signal beam. However, the choice of either beam 
might be dictated by wavelength considerations. 

Consider the situation depicted in Fig. 0. The reverse signal system hl{x,xi) = hs{xi,x) 
is regarded as a cascade of two linear systems of impulse response functions hi and /i2 with 
the object t sandwiched in between. The biphoton amplitude can thus be written as 

i.ixi,x,)=Jdx'tix')hixi,x')h,ix,,x% (5) 

where ^3 is the impulse response function of a system composed of a cascade of the reverse 
of system /i2, an aperture Ep, and the system hi, and is given by 

h,(.„.-) = ld.EAx)h.Mh,(.;.). (6) 

Equation @ states that the overall system is composed of an illumination system hi illumi- 
nating the object t, followed by an imaging system h^, which is dependent on /i2, hi, and Ep, 
in accordance with Eq. @. The uniformity of the illumination system hi and the resolution 
of the system /13 determine the quality of the overall imaging system. 

In this configuration, then, the conditional coincidence rate, obtained using Eqs. (P, 
(I), and d), is 

4'\^2) = I / dx't{x')hi{0,x')h3{x2,x')\'. (7) 

This system is mathematically equivalent to a coherent optical system where the object is 
modulated by the illumination distribution hi{0,x') and transformed by a linear system of 
point spread function h^{x2,x'), followed by a squarer, viz. the bilinear transformation of 
Eq. 

On the other hand, using Eqs. (P, (|^), and (^), the marginal coincidence rate, measured 
when Di collects photons from all points in its plane (i.e., acts as a bucket detector), is given 
by 

I^'^\x2)= j dxi\ j dx't{x')hi{xi,x)h^{x2,x')\^ 
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The quadratic transformation of the object t{x) in Eq. is clearly the mathematical 
equivalent of the bilinear transformation in Eq. (^5), representing a partially coherent 
imaging system. The function g{x', x") is given by 

g{x',x") = J dxihl{xi, x')hi{xi, x"), (9) 

and plays the role of the coherence function of the field. 

Two limiting forms of g{x',x") are of interest. The first is g{x',x") = 6{x' — x"), which 
leads to 

= / dx'\tix')\^\h{x2,x')\\ (10) 

which is the equivalent of an incoherent imaging system [Eq. ( [A3| )]. The other limit is 
g{x',x") = f*{x')f{x"), in which case 

I(^)(X2) = I / dx't{x')fix')hix2,x')\'', (11) 



which is the equivalent of a coherent imaging system [Eq. (^3)]. We can achieve the first 
limit by using a 2 — / (Fourier transform) system or a 4 — / (imaging) system for hi followed 
by a bucket detector. Moving the bucket detector in the 2 — / system away from the back 
focal plane or changing the area of the detector would lead to a gradual transition from 



coherent to incoherent imaging. This was performed experimentally in Ref. [ ||2T|] where 
the change from coherent to incoherent imaging, achieved by changing the detector aperture 
size in one beam, was observed by monitoring the loss of the fringe visibility of a double-slit 
placed in the other beam. 

We now examine a few examples based on this configuration that manifest its usefulness. 



1. Example: Fourier transform imaging 

Consider the system illustrated in Fig. |^. The signal and idler systems are unfolded for 
sake of clarity. We assume for simplicity in this and in the following examples, except if 
indicated otherwise, that the pump and the NLC are of infinite transverse extent. The signal 
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arm includes the object, t, and the system /i2 is nothing but free space propagation at the 
signal wavelength for a distance dg. The idler system comprises a lens of focal length / 
at a distance di from the NLC, and a distance ^2 from the detection plane, as shown in Fig. 
^. Free-space propagation for the distance dg at As, and for the distance di at Aj, may be 
substituted by free-space propagation for an equivalent distance di at wavelength A,, where 

di = di + ds^. (12) 

If we take d2 to be equal to the focal length of the lens /, and also choose di and dg such 



that di = f according to Eq. ([T4 ), then the system becomes a Fourier transform system 
with impulse response function 

h{x2, x') = exp(j27r[^ + ^ ])exp{~j^X2x'). (13) 

If we now take the illumination system hi to be uniform, so that hi{0,x') = 1, then 
the overall system is a Fourier transform system when the conditional coincidence rate is 
considered. Equations ([^) and (|TB|) then yield 



li>\^2) = \n^X2)\\ (14) 

where T{q) is the Fourier transform of t{x). The system simply generates the diffraction 
pattern of the object distribution. 

2. Example: Ideal single-lens imaging 

In the same configuration depicted in Fig. |^, we may choose the distance di, calculated 
according to Eq. (p!2D, and the distance ^2 to satisfy 

which is the geometrical-optics imaging equation of a thin lens of focal length /. In this 
case the impulse response function of the system /13 becomes 
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h{x„x') = exp(j27r[f + ^^])exp(j^[l - h)6{x, - Mx'), (16) 

As Aj Aja2 M 

where M = — ^ is the magnification of the imaging system. If the illumination system hi 

(2) 

is uniform, then the conditional coincidence rate /g {X2) is proportional to the magnified 
object intensity transmittance, |t(ff)p. The marginal coincidence rate I^'^\x2), however, is 
I^'^\x2) = fi'lffj If )P? where g{x',x") is given by Eq. (^. If g{x',x') is uniform over 

an area larger than that of the image, I^'^\x2) becomes proportional to the magnified object 

(2) 

intensity transmittance as is the case for Iq (3:2). 

Note that the lens may equivalently be put in the signal beam and the distances read- 
justed so as to satisfy a condition similar to Eq. (^). The system developed by Pittman et 
al. is an example of this case in which the object is placed directly in the plane of Di 



so that hi{xi,x') = 5{xi — x'). 

3. Example: Lens in the pump 

We now study another example where we manipulate the pump beam and place the 
object in either the signal or idler paths. An example of this configuration is the system 
proposed by Belinsky and Klyshko and demonstrated experimentally by Pittman et al 



The configuration is shown schematically in Fig. ^ A plane-wave pump beam is focused 
by a lens of focal length / , located at a distance d < f from the NLC The pump wave front 
now has a radius of curvature R = f — d aX the NLC entrance, and consequently acts as a 
lens or spherical mirror in the advanced wave interpretation. The signal system is comprised 
of free space propagation for a distance di at As, followed by the object t and an optical 
system hi. The idler system is simply free space propagation for a distance ^2 at Aj. 

If the following relationship is satisfied by the various distances and wavelengths: 

Xsdi Xid2 XpR 

then /i3(x2,x'), from Eq. (P), is 
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h{x„x') = exp(j27r[^ + |])exp(j^[l - j^Mx^ - Mx'), (18) 

This is the impulse response function of an imaging system of magnification M = 
and Eq. ([T7|) is the corresponding imaging equation. Note the similarities between Eqs. 
(p!5|) and (p!7[) , despite the fact that the lens is in the signal beam for the former and in the 
pump beam for the latter. 

In the degenerate case where the signal and idler frequencies are equal, Xg = Xi = 2Xp, 
this imaging equation (Eq. (0)), becomes 

112 

This is the imaging equation of a spherical mirror of radius of curvature R, or is the 
geometrical-optics imaging equation of a thin lens of focal length ^. Both l'^\x2) and 
-^'•^''(^2) are proportional to the magnified object intensity transmittance |t(||)P if hi is 
uniform. 



B. Object in Both Signal and Idler Beams 

If the signal and idler systems are identical and the object t is placed at the same 
location in each, as shown in Fig. ^ then we may substitute hs{xi,x) = hi{xi,x) = 
J dx'hi{xi, x')t{x')h2{x' , x) in Eq. to obtain 

iIj{xi,X2) = J J dx'dx"t{x')t{x")ijjc{x',x")h2{xi,x')h2ix2,x"), (20) 

where 

ijjcix' ,x") = J dxEp{x)hi{x' , x)hi{x" , x). (21) 

Comparing Eq. ( pO]) with Eq. ( |Al| ) shows that ip{xi, X2) is analogous to a partially coherent 
imaging system, where ipdx'jx") plays the role of the correlation function of the field. In 
this case, though, in accordance with Eq. ([I|) the biphoton rate G^'^\xi, X2) is a fourth-order 
nonlinear transformation of t. 
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This system has been used to test the complementarity of coherence and entan- 
glement with the change of transverse size of the pump beam, where t was taken to be a 
double-slit. From Eqs. ( pOD and (|21| ) it is clear that for a small source the biphoton ampli- 
tude factorizes into a function of xi and another function of X2 (coherence), while it is not 
factorizable (i.e., entangled) for a large pump-beam size (entanglement) [BS]. 



C. Object in the Pump Beam 

In another imaging configuration, the object is placed in the pump beam as illustrated 
in Fig. |[ Equations (|I]) and (|) give 

G^^\xuX2) = I J dxt{x)hs{xux)hi{x2,x)f (22) 

provided that Ep{x) is uniform over the object. Many possibilities for imaging based on 
Eq. ( p^ ) can be envisioned. For example, if both hg and hi are 2 — / systems, the result is 
proportional to the squared magnitude of the Fourier transform of t. In another example, if 
hs{0, x) = 1, then 

Io\x2) = I / dxt{x)hi{x2,x)\\ (23) 

and the behavior is that of a coherent imaging system. The object is illuminated at the pump 
wavelength, while the measurement is made at the much longer signal and idler wavelengths. 

In a third example in which the signal and idler systems are identical, and the coincidence 
is measured at the same position, by use of a detector sensitive to the arrival of photon pairs 
(a two-photon absorber, for example), then 

G^^^\xi,xi) = I J dxt{x)hl{xi,x)\'^. (24) 

Again, the mathematical structure is that of a coherent imaging system. 

An interesting modification to this configuration would be to add a 2 — / system between 
the object and the crystal. In this case, Ep{x) = T(27ry^), where fo is the focal length of 
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the lens before the crystal. If the object is not symmetric in then its Fourier transform is 
a complex function. Yet the phase distribution of the object's spatial spectrum is not lost, 
since the three-wave interaction process in the NLC is coherent. If, in addition, we take 
both the signal and idler configurations to be 2 — / systems, the biphoton rate becomes 

G^'\x„x,) = \t{[^ + {^]fA,)\'. (25) 

^sjs ^iji 

where fs and fi are the focal lengths of the 2 — / systems in the signal and idler beams, 
respectively. In this case, in accordance with Eq. (|), Iq'^ provides a magnified image of the 
object t with a magnification factor of 



D. Object is the Detector 

In yet another imaging modality, illustrated in Fig. |^, the object is a two-photon ab- 
sorber; it is thus a detector with quantum efficiency proportional to its absorbance ^2- The 
biphoton rate, in this case G^'^\xi,Xi), is registered by some response of the object, such as 
emitted photoelectrons or fluorescence |^3|. The signal/idler optical system may, for exam- 



ple, be a single-lens imaging system or a scanning system, as in scanning confocal microscopy 
g. From Eqs. (|l]) and (D 



G^'\xux^)=h{x,)S{x,), (26) 

where S{xi) is an object illumination function. One choice for 5' would be a very narrow 
function , which would serve to sample the object in the transverse plane. This could be 
achieved, for example, by taking a pump of large transverse width and 2 — / signal and idler 
systems. The size of the pump and the aperture of the lens limit the transverse resolution. 

There are other sensible choices for the illumination function S. These can be imple- 
mented through either the pump profile, or the system impulse response function, or both. 
The object t2 would then be extracted by dividing the observed biphoton rate by S. The 
object can also be uniformly illuminated by using a large pump and a 4 — / imaging system, 
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in which case S becomes almost constant over a large portion of the object. We have studied 
this system elsewhere and compared the longitudinal and transverse resolutions to those of 
other schemes of microscopy that utilize classical light |f45|| . 

E. Objects in Signal, Idler, and Pump Beams: Image Triple Correlation 

Because the biphoton optical system is based on three-wave mixing, it inherently depends 
on three image distributions and therefore offers a number of unique options for optical-image 
processing. For example, if 4 — / systems with aperture functions tg and ti are placed in the 
signal and idler beams, respectively, and a third object is placed in the pump beam such 
that the field at the entrance to the NLC is tp then Eqs. (0) and @ yield 

dxtp{x)T,{--{x - x,))T,i--{x - (27) 

■^sjs \Ji 

where Ts{q) and Ti{q) are the Fourier transforms of and and fs and are the 

focal lengths of the signal and idler 4 — / systems, respectively. This equation represents the 
magnitude of the triple correlation of the three functions tp, T^, and Tj. Triple correlation 
is useful in a number of signal-processing applications. Of course, if one of these three 
functions is uniform, the operation becomes ordinary correlation. 

One application for this configuration could be system identification and coded-aperture 
imaging. In this application, a linear, shift-invariant optical system is to be identified, i.e., 
its impulse response function measured. The system may be placed in one of the two down- 
converted beam (say the signal) while a set of known reference systems are placed, one 
at a time, in the other beam (the idler) as the coincidence rate is measured. The set of idler 
systems are also assumed to be linear, shift-invariant with impulse response functions 

hi{x2,x) = hnix - X2), n = l,2,...iV. (28) 

Such systems may be generated by the use of a bank of apertures (filters). Since the unknown 
system is shift-invariant, its impulse response function is hs{xi,x) = hs{xi — x), so that by 
virtue of Eqs. (|l|) and (^ the biphoton rate measured at xi = X2 = is given by 
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Cn = G'(^)(0,0) = I j dxhs{-x)hn{x)\\ (29) 

assuming the pump distribution to be uniform. 

If {hn{x)} form a complete set of orthonormal functions, then the measured coefficients 
{Cn} are simply the squared magnitudes of the coefficients of an expansion of the unknown 
function hs{—x) in this basis. Under special conditions, the phases can be retrieved, and 
the function hs{x) completely reconstructed 

In the special case for which hn{x) = 6{x + Xn), so that the idler field is sampled at 
positions x„, Eq. ( ^Of ) yields 

Cn=\hs{Xn)\\ (30) 

in which case the measured coincidence rates provide samples of the magnitude of the impulse 
response function. A scanning system can therefore be used to determine |/is(x)|. 



III. RESOLUTION OF BIPHOTON IMAGING 

In all of the configurations studied in the previous section we assumed a thin NLC and 
a narrow biphoton spectral bandwidth. Under these assumptions the imaging resolution of 
all the configurations is determined by the apertures placed in the system (including those 
placed in the pump beam). When these apertures are not accounted for, we obtain results 
reminiscent of classical geometric optics, such as the imaging formulas in Eqs. (|l^) and (|17|). 



These geometric-optics results are typical of the work that has been carried out to date in 



entangled-photon imaging ||31| , |32| . 

One of the advantages of our formalism is to facilitate deriving the analog of wave-optics 
results for such systems when all the physical parameters of the optical arrangement are 
accounted for, using straightforward calculations similar to those of classical wave optics 
PJ33]. In this section we examine the effect of the various parameters of entangled-photon 



imaging systems on the imaging resolution. 
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We take the width of the image formed by a point object t{x) = 5{x) in the marginal 
coincidence rate [Eq. (|^)] as a measure of the resolution of the entangled-photon imaging 
system. Another definition of resolution may be based on the conditional coincidence rate 
[Eq. i)]. 

We begin by modifying the principal imaging equations [Eqs. (|I|) and (^] by taking into 
consideration the thickness of the NLC and the biphoton spectral bandwidth. We assume 
a monochromatic pump beam of angular frequency Up and transverse distribution Ep{x) 
at the entrance to a NLC of thickness i. The coincidence rate G^'^\xi,ti; X2,t2), with the 
detection times of Di and D2 now made explicit, is given by 

G^^\xi,ti;x2,t2) = \'i/j{xi,ti;x2,t2)\^. (31) 

Here i/j{xi, ti, X2, ^2) may be written in terms of a biphoton spectral amplitude ip{xi, X2; ujs) 
via 



'ilj{xi,ti;x2,t2) = exp{-iujpti) J^dusexp{-ius{ti - t2))^(a;i, 0:2; a;^), (32) 

where Q is the biphoton spectral bandwidth and 'i/j{xi,X2;uJs) is given by @] 

ijj{xi,X2;iOs) = J J dqsdqiA{qs,qi;uJs)Hs{xi,qs;iOs)Hi{x2,qi;ujp - iOs), (33) 

and the dispersion of the optical systems has been made explicit in the signal and idler 
transfer functions Hg and Hi, which are Fourier transforms of the impulse response functions 
kg and hi (with respect to the second argument), respectively. The quantity A{qs,qi]ujs) in 
Eq. (|33D is given by 

A{qs, qf, uJs) = Ep{qs + qi)i{qs, qu ^s); (34) 

here q is proportional to the transverse component of the momentum vector (it is the spatial 
frequency in the transverse plane), Ep{q) is the Fourier transform of Ep{x), and ^{qs, qf, ujg) 
is a phase-matching function given by 

£ i 

({q.,, qf, Us) = £sinc(— Ar)exp(-j-Ar); (35) 
15 



and Ar = rp(g^ + Qi, ujp) - rs{qs,uJs) - ri{qi, iOp - Us)] rj{q, u) = Jn]^ - q^, j = p, s, and i, 
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where nj is the NLC index of refraction for the polarization and frequency of the j field. 

In most cases the detectors may be considered to be slow (i.e., their response time is 
large with respect to the inverse of the bandwidth of the system, which is a reasonable 
assumption for available photodetectors), and thus they measure a coincidence rate that is 
averaged over a long time interval. The resulting time averaged coincidence rate is |53I 



(36) 

showing that the time averaged coincidence rate for a slow detector is an incoherent sum of 
the biphoton spectral amplitudes over the bandwidth of the system. The spectrum of the 
down-converted biphotons can be quite large and the dispersion of the optical components 
must be considered carefully just as dispersion must be in ultrafast pulsed optics. 
We also define conditional and marginal time-averaged coincidence rates as 

C{x2) = J dxiC{xi,X2), (37) 

C7o(x2) = C(0,X2), (38) 

respectively. It is obvious that when only a narrow spectral bandwidth is considered, C{x2) 
and Co{x2) coincide with /(X2) and Io{x2), respectively. 

We now proceed to study the resolution of a representative configuration considered in 
Section 2: object in the signal (or idler) beam. The biphoton spectral amplitude of this 
system, illustrated in Fig. ^ now taking into consideration the thickness of the crystal and 
spectral bandwidth of the system, is given by 



ip{xi,X2]UJs) = j dx't{x')hi{xi,x';Us)h3{x2,x';Us), (39) 

where 

.,(.,.;.,^//.,..„A„.„.;.,.,.,.;.,H,.„.;.,-.,. (40, 
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We assume throughout that the object is thin and non-dispersive. To determine the res- 
olution of this imaging configuration we take t{x) = S{x), whereupon Eq.(^) becomes 

■iIj{xi,X2;uJs) = hi{xi,0;Us)h3{x2,0;uJs), and consequently 

C{x2) = / duJs\h3{x2,0;Us)\'^g{uJs), (41) 
Jn 

Co{x2) = / dus\h{x2,0;uJs)\'^go{uJs), (42) 

where g{uj) = J dx\hi{x,0;LL!)\'^ and go{uj) = |/ii(0, 0; Note that the system hi affects 
the imaging resolution only through introducing an effective spectral bandwidth that may 
be ignored if it is larger than that of h^. 

As a concrete example we consider the system examined in Section 2. A. 2, which is the 
second example of object in the signal beam configuration, namely ideal single-lens imaging, 
illustrated in Fig. ^. We assume, at first, a plane wave pump, so that /i3(x2, 0; Ug) simplifies 
to 

hs{x2, 0; Us) = J dqsiiQs, -Qs] uJs)H2{0, Qs] ujs)Hi{x2, -Qs] ujp - ujs). (43) 
In this example, the transfer functions of the systems h2 and hi are given by 

H2{0, Qs] Us) = exp{jksds)exp{-j^^), (44) 

Jxi'OC^ d'(^ Ix:'3j'2 

Hi{x2, -Qs] uJp - Us) = exp{jki{di + d2))exp{j^--^)exp{-j^)Pg{qs + ^— ), (45) 

2rf2 2ki d2 

where Pg{q) is the Fourier transform of p(a;)exp(j^Y^[^ — j]) with respect to x, and p{x) is 
the lens aperture. Substituting Eqs. (0), (@) and (^51) into Eq. (|^) we obtain /i3(x2, 0; ujs) 
which we then use in Eqs. ( PD and (^) to estimate the resolution. 

There are two techniques to implement this system in an actual setup. In one technique 
the NLC is adjusted for non-collinear SPDC, and one beam (usually chosen by a pinhole) 
is directed into the system hg and the other into hi. Another technique is to adjust the 
NLC for collinear SPDC and then separate the two photons comprising the biphoton. In 
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type-II SPDC (where the signal and idler photons have orthogonal polarizations) one can 
use a polarizing beam splitter to separate the biphoton. On the other hand, in type-I SPDC 
(where the signal and idler photons have the same polarization) the use of a non-polarizing 
beam-sphtter will separate the pair into the two output ports of the beam splitter in 50% 
of the trials, and send the pair together into one output port in the rest of the trials. In 
the latter case, the trials do not contribute to the coincidence measurements carried out by 
the detectors Di and D2 together with the coincidence detection circuit, and thus may be 
ignored. 

Assuming a thin NLC, narrow spectral bandwidth, a plane-wave pump, and degener- 
ate collinear down-conversion (where the signal and idler photons are separated using the 
method outlined above), one obtains the familiar diffraction pattern of a diffraction-limited 
imaging system. For a rectangular lens aperture of width D and focal length / the result 
is C(x2) oc |sinc(2]f^)P, where = is the the F-number of the lens and Aq = 2Ap is 
the wavelength of the degenerate down-converted photons. This is the best one can obtain; 
we demonstrate in the following that relaxing any of the restrictions indicated above will 
degrade the resolution. 

Our calculations have been carried out using a /3-barium borate (BBO) NLC that is 
illuminated with a pump of wavelength Ap = 325 nm (which corresponds to the ultraviolet 
line of a He-Cd laser), with a cut-angle of 36.44° that corresponds to degenerate collinear 
type-I SPDC. Increasing the cut-angle beyond 36.44° yields non-coUinear degenerate SPDC 



while decreasing the cut-angle below this value yields collinear non-degenerate SPDC ||2^ . 

We first consider the effect of the finite thickness of the NLC. One effect is that the 
distance di, used in the imaging formula presented in Eq. (p!5|), is modified to become 



di = di + + ieq, (46) 

in contrast to that given in Eq. (0). The quantity leq is an equivalent length for the NLC 
that is related to the physical length I by 

4, = ^(^ + ^). (47) 

2 A,- Uj 



For the degenerate case (A^ = Aj = Ao) this expression simphfies to 4g = ^7-, where Uo 
is the index of refraction of the NLC at the degenerate wavelength. In other words, the 
thickness of the NLC must be accounted for in calculating the distances in the experimental 
arrangement in order to satisfy the imaging formula of Eq. (|15]). 

Figure ^ shows C{x2) for several values of i, assuming degenerate collinear down- 
conversion and narrow spectral bandwidth. The distances in the configuration are chosen 
such that di = d2 = 2/ (taking into account the effect of £ on di according to Eq. (^Bf)), 
which corresponds to an imaging system of unity magnification. For £ = 0.1 mm one obtains 
the diffraction limited distribution corresponding to the thin NLC case. When £ increases 
the distribution widens, signifying a loss of imaging resolution, as is evident for the £ = 1- 
mm and £ = 10-mm curves. This result may be easily understood when one considers the 
fact that the NLC acts as a spatial filter through the phase-matching function ^{qs,qi]UJs) 
defined in Eq. (^). The collinear SPDC case corresponds to a low-pass spatial filter with 
a cut-off frequency that is inversely proportional to £ and hence the resolution degrades as 
the NLC thickness increases. 

The spectral bandwidth of the system has a similar effect on the imaging resolution. 



which decreases with increased bandwidth. Figure |T0] shows C{x2) for several values of 
p = These plots were obtained for a NLC of thickness £ = 1 mm, collinear SPDC, and 
a plane-wave pump. According to Eq. ( PB| ) di is a function of wavelength (and so is £eg via 
Eq. d^)), so that only one pair of signal/idler wavelengths satisfy the imaging formula in 
Eq. (plsl) . All biphotons with other signal/idler pair wavelengths are defocused, and hence 
their contribution to C{x2) leads to a reduction in resolution. The plots in Fig. |l^ were 
obtained assuming that Eq. (|l^) is satisfied by the degenerate signal/idler wavelengths, and 
that at these wavelengths di = d2 = 2f. 

Finally, the finite transverse width of the pump field also degrades the resolution. This 
can be understood by noting that smaller pump size reduces entanglement of the signal and 
idler photons [p3]. As a result, the quantum state of the light emitted by the NLC becomes 
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separable and thus C{x2) and Co{x2) simply become the intensity of the idler beam (which 



depends on h^), but are independent of the signal beam No information about the 

system hg, which includes the object to be imaged in this case, may be extracted from the 
measurements carried out in the idler beam. 

Figure O shows plots of C{x2) for various values of the transverse width of the pump, 
denoted B. The calculations were performed taking £ = 1 mm, assuming collinear degenerate 
SPDC, and the presence of narrowband spectral filters in the system. Distances were chosen 
such that di = d2 = 2f. 



IV. CONCLUSION 

We have presented a Fourier-optics analysis of various imaging configurations using the 
unique features of spontaneous parametric down-conversion (SPDC) as a two-photon source. 
SPDC is a three-wave mixing process; the pump, signal, and idler are coupled through the 
phase-matching conditions. We investigated several imaging and image-processing configu- 
rations that utilize the quantum correlations among these three fields. Our formalism was 
also used to study the resolution of these entangled-photon imaging configurations. 



APPENDIX: THE OPTICAL BILINEAR TRANSFORMATION 

We present a brief overview of the theory of classical imaging in the framework of the 
bilinear optical transformation. The equations are formulated in such a way so as to facilitate 
convenient comparisons with the two-photon and biphoton cases presented in the text. 

Because of the quadratic relation between the optical field and the optical intensity, 
imaging systems are typically described by a bilinear transformation [^. A general bilinear 
transformation is expressed as 

g{xi)=J J dx'dx"f*{x')f{x")q{xi;x',x"), (Al) 
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where f{x) is the input function, q{xi; x', x") is the double impulse response function (DIR), 
and g{x) is the output function. In general f{x) is complex, but g{x) is guaranteed to be real 
when the symmetry condition q{xi\x'^x") = q*{xi] x" , x') is satisfied. The DIR completely 
characterizes the bilinear system. This transformation represents, in general, the imaging 
system depicted in Fig. O. The input function f{x) represents the transparency t{x); the 



DIR is a combination of the second-order correlation function of the illumination G^^\x', x") 
and the impulse response function h{xi, x) of the linear optical system; and the output g{x) 
represents the intensity measured by the optical detector. 

In the ideal case q{xi; x', x") = 6{xi — x')6{xi — x"), whereupon g{xi) = so that 

the system is a squarer with zero spread. When the DIR factorizes in the form q{xi; x\ x") = 
h*{xi,x')h{xi,x"), the output of the system is given by 



g{xi) = I J dx'f{x')h{xi,x')\\ (A2) 

Equation ( |A2| ) is easily recognized as the output intensity of a coherent imaging system with 
impulse response function h{xi,x') and an input complex field f{x'). When the DIR takes 
the form q{xi]x' ,x") = h*{xi,x')h{xi,x")6{x' — x"), we obtain 

g{xi)= j dxy{x')\^\h{x,,x')\\ (A3) 

which is the output of an incoherent system with point spread function x')\'^ and input 

intensity \ f{x')\^. In general, partially coherent imaging can be represented by a bilinear sys- 
tem with a DIR given by q{xi\ x\ x") = 7(2;', x")h*{xi, x')h{xi, x"), where 7(2;', x") represents 
the correlation function of the input light, and h{xi,x') is the coherent impulse response func- 
tion. When 7(x', x") = 1, we recover coherent imaging whereas when j{x', x") = 6{x' — x"), 
we recover incoherent imaging. 

Entangled-photon imaging, like partially coherent imaging, is described by a bilinear 
system, with partial entanglement assuming the role of partial coherence |^ . 
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FIGURES 

FIG. 1. Biphoton imaging using photon pairs generated by spontaneous parametric 
down-conversion. NLC stands for nonlinear crystal; Di and D2 are single-photon detectors at lo- 
cations xi and X2, respectively; G^'^^xi, X2) is the biphoton rate; hp{x,x'), hs{xi,x), and hi{x2,x) 
are the impulse response functions of the optical systems placed in the paths of the pump, signal, 
and idler beams, respectively. 

FIG. 2. Object in the signal beam configuration. Ep is the pump field at the entrance to the 
NLC; hi{xi,x') and h2{x',x) are the impulse response functions of the optical systems placed in 
the signal beam; hi{x2, x) is the impulse response function of the optical system placed in the idler 
beam; t{x') is the object to be imaged, placed in the signal beam. 

FIG. 3. Object in the signal beam configuration of examples 1 and 2 displayed in an unfolded 
picture. Ep, hi{xi,x'), and t{x') are the same as in Fig. |2|; / is the focal length of a lens placed in 
the idler beam. See text for details. 

FIG. 4. Object in the signal beam configuration of example 3 displayed in an unfolded mode. 
hi{xi,x'), and t{x') are the same as in Fig. |2|; a lens is placed in the pump beam and is represented 
here by the dotted lens of focal length /. See text for details. 

FIG. 5. Object in both signal and idler beams configuration. hi{xi,x'), h2{x', x) are the impulse 
response functions of the optical systems placed in the path of both the signal and idler beams; 
and t{x') is the object to be imaged. 

FIG. 6. Object in the pump beam configuration. hs{xi,x), hi{x2,x) are as in Fig. ||; and t{x) 
is the object to be imaged placed in the pump beam . 

FIG. 7. Object is the detector configuration. hi(xi,x), is the impulse response of the optical 
system placed in the signal and idler paths; and D is a two-photon detector at location xi. 
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FIG. 8. Configuration for triple correlation. h2s{x',x) and his{xi,x') are the impulse response 
functions of the optical systems placed in the signal beam; h2i{x", x) and hii{x2, x") are the impulse 
response functions of the optical systems placed in the idler beam; tp{x), ts{x'), and ti{x") are the 
three objects to be correlated. 

FIG. 9. Effect of NLC thickness £ on the imaging resolution of object in the signal beam con- 
figuration. Plots of normalized time-averaged marginal coincidence rate C{x2) versus detector D2 
location X2 normalized with respect to Xc = 2XoF^, for i = 0.1, 1, and 10 mm; Aq = 650 nm, 
F# = 5. 

FIG. 10. Effect of the biphoton bandwidth fl on the imaging resolution of object in the signal 
beam configuration. Plots of normalized C{x2) versus X2 normalized with respect to Xc = 2Ao-P# 
are shown for /9 = ^ = 0.001, 0.01, and 0.02; = 650 nm, F# = 5. The NLC is adjusted for 
collinear SPDC and is of thickness 1 mm. 

FIG. 11. Effect of the transverse width of the pump B on the imaging resolution of object in 
the signal beam configuration. Plots of normalized C{x2) versus X2 normalized with respect to 
Xc = ^XoF^ are shown for B = 2,1, 0.5, and 0.1 mm; Aq = 650 nm, = 5. 

FIG. 12. Classical partially-coherent imaging. The quantity {x',x") is the second-order 
correlation function of the optical field; t is the object to be imaged; h{xi,x) is the impulse 
response function of the imaging system; D is a detector placed at position xi that records the 
intensity I(xi). 
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